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Abstract 

Combining data on unpolarized and polarized inclusive proton structure functions, we perform 
the first detailed study of quark-hadron duality in individual helicity- 1/2 and 3/2 virtual photo- 
production cross sections. We find that duality is realized more clearly in the helicity-1/2 channel, 
with duality violating corrections < 10% over the entire nucleon resonance region, while larger, 
< 20% corrections are found in the helicity-3/2 sector. The results are in general agreement with 
quark model expectations, and suggest that data above the A resonance region may be used to 
constrain both spin-averaged and spin-dependent parton distributions. 



I. INTRODUCTION 



The duality between quark and hadron descriptions of physical observables reveals a 
fascinating connection between the physics of quark confinement at low momentum scales 
and asymptotic freedom at large momenta. The striking manifestation of this duality in 
inclusive electron-nucleon scattering, first observed [l| even before the advent of QCD, has 
motivated considerable effort in recent years to explore this phenomenon empirically, as well 
as to understand its origins theoretically (for a review see Ref. j^J). 

From a practical perspective, the quantitative demonstration of the similarity between the 
structure functions measured in the nucleon resonance region and those in the deep inelastic 
continuum, at higher energies, opens up the intriguing possibility of using resonance region 
data to provide constraints on leading twist parton distribution functions (PDFs). Attempts 
to utilize this connection have begun to be explored in recent global PDF fits js| (see also 
Ref. {4 1), where data on the unpolarized proton and deuteron F 2 structure functions at final 
state hadronic masses W as low as ~ 1.7 GeV have been used to extend determinations of 
PDFs to larger values of the Bjorken scaling variable x = Q 2 /2Mu, where Q 2 and v are the 
four-momentum squared and energy transferred to the proton, and M is the proton mass. 

The availability of high-luminosity electron beams at Jefferson Lab has enabled high- 
precision measurements of various structure functions to be made over the past decade. 
These data have now firmly established the existence of duality in the proton F 2 and Fl 
structure functions Q"8L and have provided tantalizing glimpses of its s pin and flavor depen- 
dence in polarized 9|, [l0| and semi-inclusive scattering measurements llfl. Recently a new 



method 



12j was used to extract also the neutron F 2 structure function from inclusive proton 



and deuterium data in the nucleon resonance region [13|], leading to the first quantitative 
determination of duality in the neutron's unpolarized structure functions. This observation 
suggested that duality is indeed a general feature of the resonance-scaling transition, and 
not due to accidental cancellations of quark charges [ljj]. 

Because of the considerably larger data base of spin-averaged cross sections than polar- 
ization asymmetries, duality in the spin-dependent g\ and g 2 structure functions has not 
yet been established to the same precision as for the unpolarized Fx s 2 structure functions. 
An additional complication arises from the fact that for spin-dependent quantities one deals 
with differences of cross sections, which are not restricted to be positive. For example, in 
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the A resonance region the g\ structure function of the proton, especially at low Q 2 , is large 
and negative, while for the same x and higher Q 2 (hence higher W 2 = M 2 + Q 2 (1 — x)/x) the 
structure function measured in the deep inelastic region is positive. Such strong violation 
of duality would limit the use of resonance region data to constrain spin-dependent PDFs. 
Furthermore, the neutron gi structure function changes sign as a function of x, so cannot 
be used to study ratios of resonance to deep inelastic structure functions or the relative size 
of duality violations. 

On the other hand, duality violation may be less severe, even at low W, in individual 
virtual photoabsorption helicity cross sections, defined by projecting the total spin of the 
7*-proton center of mass system along the photon direction. The helicity-1/2 projection, 
CT1/2, represents the cross section for equal initial and excited hadronic state helicities, while 
the helicity-3/2 projection, (X3/2, involves a change of the hadron helicities by two units. The 
sums and differences of the helicity cross sections, which are positive definite, correspond to 
the unpolarized and spin polarized structure functions, respectively. 

At high Q 2 and W 2 the helicity cross sections are proportional to the positive and negative 
helicity PDFs, q ± (x,Q 2 ), which describe the distribution of quarks with spin parallel or 
antiparallel to that of the nucleon. As with the cross sections, the helicity PDFs are defined 
to be positive, and in a way represent more fundamental objects than the spin-averaged and 
spin- dependent PDFs. In fact, from perturbative QCD arguments one can make definite 
oredictions for the behavior of the helicity PDFs in the limit x — > 1, with q~ /q + ~ (1 — x) 2 
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161 ] . These predictions can be tested by studying the asymptotic x dependence of helicity 



cross sections, which is difficult, however, because of the rapidly decreasing rates as x — > 1. 
Indeed, since large x generally corresponds to low W, determining the large-x behavior of 
inclusive cross sections at any finite Q 2 will necessarily involve the resonance region. 

In this paper we perform the first detailed study of quark-hadron duality in 7*p helicity 
cross sections, by combining previously measured sets of data on inclusive spin-averaged cross 
sections and double polarization asymmetries from Jefferson Lab and elsewhere. These data 
are used to quantify the degree of duality violation in each of the three prominent nucleon 
resonance regions, as well as over the entire range W < 2 GeV. 

In Sec. [Til we begin by defining the relevant cross sections and distributions used in this 
analysis. Section II III outlines the data analysis, describing the construction of the helicity 
cross sections form separate measurements of spin-averaged and spin-dependent structure 
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functions. Results of the analysis are presented in Sec. HVt for the x dependence of helicity 
structure functions in several fixed-Q 2 bins ranging from Q 2 = 1.7 to 5 GeV 2 , as well as 
for integrals over the various resonance regions. Comparison of the resonance data with 
parametrizations of data at higher energies then allow the first determination of the extent 
to which duality holds in helicity cross sections over this range. We also compare our 
findings with quark models that predict specific patterns of duality violation in structure 
functions. Finally, in Sec. |V] we draw some conclusions from this analysis and outline its 
broader implications for our understanding of quark-hadron duality as well as its practical 
exploitation. 



II. HELICITY STRUCTURE FUNCTIONS 



The differential cross section for the inclusive scattering of a longitudinally polarized 
electron with helicity h = ±1 from a proton with polarization P z along the virtual photon 
direction can be written as 



da 



dVtdE' 



T (a T + ea L + hP z y/\- e 2 a' TT ) , (1) 



where E' is the scattered electron energy, T is the flux of virtual photons and e is the 
t_e p ll0ton pel— 0. The p—rption cross sections for transversely 
polarized virtual photons are related to the helicity cross sections ai/ 2 and a 3 / 2 by 

a T = 2 ( ai /2 + °3/2) , (2) 
a' TT = - (0-3/2 - 0-1/2) , (3) 

while the cross section for longitudinally polarized photons is given by the longitudinal 
structure function. The cross section ai/ 2 (03/2) corresponds to the spins of the virtual 
photon and proton antialigned (aligned) in the center of mass system, so that the helicity of 
the excited nucleon state after absorbing a photon is +1/2 (+3/2). Whereas the ai/ 2 cross 
section conserves the nucleon helicity, the 03/2 changes the nucleon helicity by two units. 

For convenience we define dimensionless helicity structure functions Hi/ 2 and H 3 / 2 i n 
terms of the cross sections by 

MK MK 
Hl/2 = 4^ ai/2 ' Hi/2 = 4^ a3/2 ' (4) 



where a = e 2 /47r and K = (W 2 —M 2 )/2M is associated with the choice of the virtual photon 



flux in the Hand convention [18]. The helicity structure functions can then be written in 
terms of the usual unpolarized F 1 and polarized g lt2 structure functions as 

Q 2 



H1/2 

H3/2 



F\ 



9i 
9i 



Q 2 

-^92, 



(5a) 
(5b) 



each of which is a function of two variables, typically taken to be x and Q 2 . In the limit 
where both Q 2 and W 2 are large, with x finite, the F\ and gi structure functions can be 
written, at leading order in a s , in terms of leading twist PDFs, 



q 



(6a) 
(6b) 



where q = q + + q and Aq = q + — q are the spin-averaged and spin-dependent PDFs. In 
this case the helicity structure functions become 



H1/2 
H3/2 



4 (<r + q 



(7a) 
(7b) 



so that in this limit H\i2 is determined by the q + PDFs while H3/2 is determined by the q~ 
PDFs (the antiquark distributions q + and q~ are suppressed by additional powers of (1 —x) 
compared with the quark PDFs). In the x — > 1 limit the leading behavior of the helicity 
distributions is predicted from perturbative QCD to be q + ~ (1 — x) 3 and q~ ~ (1 — x) 5 if 



the nucleon ground state wave function is dominated by its S*-wave component 
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q ~ (1 — x) 5 log 2 (l — x) if one includes orbital angular momentum 19]. 



16| . or 



III. DATA ANALYSIS 



The experimental i/1/2 and i/3/2 helicity structure functions used in this analysis were 
obtained by combining measurements of gi/F\ ratios from Jefferson Lab experiment EGlb 
(E9 1-023) in CLAS 20] with the unpolarized F\ structure function from the empirical 



5 



Christy-Bosted (CB) global fit 2l|. For the small correction from the g 2 structure func 



22]. 



tion we use the phenomenological parametrization of Ref. 

The CLAS E91-023 data set represents one of the few high-precision measurements of 
gi/Fi for the proton at moderate to large x (x > 0.15), covering a large Q 2 range from 
0.05 to 5 GeV 2 , over both the resonance and DIS regions. The empirical CB global fit uses 
measurements of inclusive inelastic electron-proton cross sections in the kinematic range of 
Q 2 < 8 GeV 2 and W between 1.1 and 3.1 GeV. The fit is constrained by high precision 
ongitudinal and transverse separated cross section measurements from Jefferson Lab Hall C 
6], unseparated Hall C measurements up to Q 2 of 7.5 GeV 2 and photoproduction data 
at Q 2 = 0. This fit was chosen because it covers a wide kinematic range and uses both trans- 
verse and longitudinal cross sections, which is particularly important for the F\ estimation. 
Due to the scarcity of g 2 measurements, especially in the resonance region, we use the phe- 
nomenological parametrization of Ref. 22j, which is developed for x > 0.02 using DIS data 
with Q 2 up to 50 GeV 2 as well as experimental results on both photo- and electroproduction 
of proton resonances. 

The statistical uncertainties for Hi/ 2 and i/3/2 were calculated from those of the g\/F\ 



measurements 



20] . The systematic uncertainties were obtained from those of gi/Fi, Fi and 



g 2 by varying these quantities within the limits given by their systematics. The resulting 
variations of i/1/2 and H 3 / 2 were then added in quadrature to obtain the total systematic 
uncertainty. Note that since H\/ 2 and H3/2 are different combinations of the same structure 
functions (differing only in relative sign, Eqs. (jSJ)), the resulting absolute uncertainties are 
the same for the two cross sections. However, since H 3 / 2 <C -ffi/2 the relative uncertainties 
will be different, with that on H 3 / 2 much greater than on Hi/ 2 . 

The theoretical H\/ 2 and H 3 / 2 structure functions were obtained by combining g\ from the 
Bliimlein-Bottcher (BB) global parametrization 23[ of spin-dependent structure functions 
with F\ constructed from the F 2 global fit of Alekhin et al. [4] and the R1998 parametrization 
of the longitudinal to transverse cross section ratio R 24j . The BB global fit 23[ is based on 
a next-to-leading order QCD analysis of the world data on polarized deep inelastic scattering, 
and includes possible higher twist contributions. The analysis finds that for both proton 
and deuteron targets the higher twist corrections to g\ are consistent with zero, within the 
large uncertainties of the data. For the g 2 structure function we therefore use the Wandzura- 
Wilczek relation 25[ with gi from the BB fit. 
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The global fit of Alekhin et al. [4J provides QCD parametrizations for both the Fx and F2 
structure functions. While the F2 fit reproduces well the available F2 data, the Fx fit shows 
some discrepancies at low W with the high-precision longitudinal and transverse separated 



cross section measurements from Jefferson Lab Hall C 



6|. We find that a good description 



of the Hall C F\ data can be obtained by using the F 2 fit from Ref. [4] with the R1998 
parametrization of R. 

IV. RESULTS 

The results for the helicity structure functions Hx/2 and H 3 / 2 are shown in Figs. Q]and|2] 
as a function of x for several fixed Q 2 values ranging from Q 2 = 1.7 GeV 2 to 5 GeV 2 . The 
data are compared with curves (labeled "theory" ) constructed from global fits to structure 
functions in the deep inelastic region at higher W, as outlined in Sec. II III above. The 
resonance region data are in excellent agreement with the global fit for the Hx/2 structure 
function for the kinematics considered. The agreement for H3/2 is also quite good overall, 
although here the "theory" curve slightly underestimates the data, especially at lower Q 2 
and in the A resonance region, where a prominent peak stands out. This can be understood 
from the fact that at low Q 2 the A contribution to Fx is positive, while that to g\ is negative, 
thereby cancelling in H 1 / 2 but reinforcing in H 3 / 2 - 

The degree to which duality holds can be quantified by considering integrals of the struc- 
ture functions over individual resonance regions, AW, 

I(AW,Q 2 ) = [ dxF(x,Q 2 ) : (8) 

J AW 



where F = Hi/ 2 , H 3 / 2 , Fi or g 1 . Following earlier data analyses {5, 8, |13] , we take for A the 
three prominent resonance regions, defined on the intervals 

• 1 st resonance region: 1.3 < W 2 < 1.9 GeV 2 

• 2 nd resonance region: 1.9 < W 2 < 2.5 GeV 2 

• 3 rd resonance region: 2.5 < W 2 < 3.1 GeV 2 

as well as the entire resonance region W 2 < 4 GeV 2 . These are shown in Figs. [3] and H] for 
the Hx/2 and Hz 12 cases, respectively, with the integrals for F\ and g\ shown in compari- 
son. Because Hy2 involves a sum of the (positive) Fx and (generally positive) gx structure 
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FIG. 1: Scaled helicity-1/2 cross section H\j2 as a function of x for various Q 2 bins. The bands 
(labeled "theory") represent a global fit to high-VF data (see text). 

functions, the H1/2 data are generally larger in magnitude than F\ and g\. Above the first 
resonance region the agreement with the global fits (shaded regions) is extremely good over 
the entire range of Q 2 considered for each of the H 1 / 2 , i*i and g\ structure functions. The 
agreement in the first resonance region is markedly worse, reflecting the strong violation 
of duality in the g\ structure function in the vicinity of the A resonance. This violation 
persists until Q 2 ~ 3 GeV 2 , above which the resonance and deep inelastic data are in better 
agreement. 

For the H^/ 2 structure function, because this involves the difference between F\ and gi, 
its magnitude is considerably smaller than that of F\. Again, duality violation is strongest 
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FIG. 2: As in Fig. [T] but for the helicity-3/2 cross section H 3 / 2 . 

in the A region, with generally good agreement between resonance and deep inelastic data 
at higher W. While the violations of duality are expected to diminish at larger Q 2 , the 
decreasing magnitude of the higher-Q 2 integrals makes it more difficult to quantify the 
violation accurately. 

To ameliorate this problem we compute the ratios of the integrals of the resonance region 
data to those of the global fits, shown in Fig. [5] for H\i 2 and H 3 / 2 . in general these ratios 
show that duality violation is stronger in the helicity-3/2 channel than in the helicity-1/2, 
with the duality violating corrections for H3/2 positive in the first resonance region and 
negative in the second resonance region. As could be expected, the differences between the 
data and theory are largest in the A region for both Hi/ 2 and H 3 / 2 . The larger uncertainties 



9 



0.04 



~i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 



▼ data 1 st 



0.05 



t — i — i — I — i — i — i — I — i — i — i — I — i — i — i — I — i — i — r 



i 1 1 1 i 1 1 1 i 
a EH gl 




GeV 



123456 123456 



Q 2 (GeV 2 ) 



FIG. 3: Integrals I (AW, Q 2 ) of the scaled helicity-1/2 structure function -Hi/2 in various resonance 
regions AW (1st, 2nd, 3rd and W < 2 GeV) versus Q 2 . For comparison the corresponding integrals 
of the gi and F\ structure functions are also shown. The bands represent a global fit to high-11^ 
data (see text). 



on the H 3 / 2 data reflects the fact that if 3/2 -C H\ii- Integrating over the entire W < 2 GeV 
region, the duality violation is < 10% for H 1/2 and < 20% for H 3/2 at Q 2 < 4 GeV 2 . This is 
considerably smaller than the corresponding duality violation found in the spin- dependent 
gi structure function 

Our results can be compared with quark model predictions for the relative strengths of 



the N — > N* transitions. In Tab. 1 these are displayed 
SU(6) P = 56+ (L = 0) and 70~ (L = 1) representations 



or if 1/2 and if 3/2 in the various 
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27j |. with each representation 
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FIG. 4: As in Fig. [3] but for the helicity-3/2 structure function H^/ 2 - 

weighted equally. The contributions from the symmetric and antisymmetric components of 
the ground state nucleon wave function enter with strengths A and p, respectively, and the 
SU(6) limit corresponds to A = p. The usual quark model assignments of the excited states 
have the nucleon and A in the quark spin-| 2 8 and quark spin-| 4 10 representations of 
56 + , respectively. For the odd parity states the 2 8 multiplet contains the states 5*11(1535) 
and Z?i3(1520), which dominate the second resonance region; the 4 8 contains the <Sii(1650), 
Di 3 (1700) and I>i 5 (1675); and the isospin-f states S 31 (1620) and D 33 (1700) belong to the 
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10 representation [281 ]. 



With the exception of the A region, the Hi/ 2 structure function is predicted to be much 
larger than the H3/2, as is borne out by the data in Figs. [1] and [2j The relatively small 
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FIG. 5: Ratio (data - theory) /theory for the scaled helicity cross sections and H 3 / 2 m various 
resonance regions AW (1st, 2nd, 3rd and W < 2 GeV) versus Q 2 . Some of the points have been 
offset for clarity. 

contribution to Hi/ 2 in the 4 10[56 + ] channel and large contribution in the 2 8[70~] channel 
suggests that the helicity-1/2 data should lie below the global fit in the first resonance region 
and above the global fit at larger W . This is generally consistent with the data in Fig. |5j 

The helicity-3/2 structure function is dominated in the resonance region by the A, with 
suppressed contributions in all other channels. Again this is consistent with the i/3/2 data 
being higher than the global fit in the A region and below the fit at larger W . The prediction 
of vanishing H3/2 for the nucleon elastic contribution reflects the dominance of magnetic 
coupling assumed in the model {28 1, which is expected to be a better approximation at high 
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TABLE I: Relative strengths of N — > N* transitions for helicity structure functions #1/2,3/2 i n the 



SU(6) quark model [26|. The coefficients A and p denote the relative strengths of the symmetric 
and antisymmetric contributions of the SU(6) ground state wave function, with the SU(6) limit 
corresponding to A = p. 
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Q 2 . 

Summing over all channels, the ratio of helicity-3/2 to 1/2 structure functions is pre- 
dicted to be H3/2/H1/2 = 2/7, which coincides exactly with the quark-parton model results 
u~/u + = 1/5 = d + /u + and d~~ /d + = 2 = d~/u~ for all x. These predictions are found to 
hold approximately at x ~ 1/3, but significant deviations are observed at larger x. Various 
scenarios for SU(6) symmetry breaking, consistent with quark-hadron duality, were consid- 



ered in Ref. 



27|, leading to specific predictions for structure function ratios in the x — > 1 



limit. The general trends of the duality violations persist even in the more realistic sym- 
metry breaking scenarios, so that the deviations from unity in Fig. [5] can be understood, 
at least qualitatively, in terms of a microscopic quark-level description. Note also that the 
quark model predictions relate to the resonant components of the data only; the presence of 
the nonresonant background washes out these predictions somewhat, especially at larger Q 2 , 
and its remarkable that the general trends of the duality violations in the various resonance 
regions nevertheless remain. 



V. CONCLUSIONS 



In this work we have performed the first detailed analysis of quark-hadron duality in 
individual 7*p helicity cross sections, utilizing recent data on inclusive unpolarized and 
polarized structure functions from Jefferson Lab. Unlike spin-dependent structure functions 
which can change sign as a function of x and Q 2 , the helicity cross sections are by definition 
constrained to be positive definite. This reduces the dramatic violations of duality seen for 
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example in the proton g\ structure function in the A resonance region, where the negative 
resonance contribution at low Q 2 makes way for a positive structure function in deep inelastic 
kinematics at large Q 2 . 

The data on the polarized and unpolarized structure functions are used to quantify the 
degree of duality violation in the helicity structure functions in each of the three prominent 
nucleon resonance regions, as well as over the entire range W < 2 GeV. We find that 
duality is realized more clearly for helicity- 1/2 structure function if 1/2 than for the helicity- 
3/2 function H3/2, with the duality violating corrections in the latter positive in the first 
resonance region and negative in the second resonance region. The duality violations are 
largest in the A region for both if 1/2 and H^/ 2 - 

Over entire resonance region, the duality violating corrections are < 10% (and negative) 
for H1/2 and < 20% (and positive) for H3/2 at Q 2 < 4 GeV 2 , which is rather smaller than 
the corresponding duality violation found in the spin-dependent gi structure function. The 
patterns of duality violation are in general agreement with expectations from quark models 



based on spin-flavor symmetry 2a, |27 |. 



Our results suggest that data above the A resonance region could be used to constrain both 
spin-averaged and spin-dependent parton distributions. This lends support to recent efforts 
to broaden the kinematic coverage in global fits of unpolarized PDFs 3] by lowering the Q 2 
and W 2 cuts, and to extending these efforts to the polarized sector. Moreover, it raises the 
interesting possibility of performing global fits of helicity PDFs q + and q~ directly, rather 
than reconstructing these from separate unpolarized and polarized PDF analyses; such an 
enterprise would demand a consistent analysis of combined cross section and polarization 
asymmetry data along the lines presented in this work. 
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